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Ab-initio study of the anomalies in the 
He atom scattering spectra of H/Mo (110) and H/W (110) 

Bernd Kohler, Paolo Ruggerone, and Matthias SchefHer 
Fritz-Haber-Institut der Max- Planck- Ges ells chaft, 
Faradayweg 4-6, D-I4 195 Berlin- Dahlem, Germany 
(February 8, 2008) 

Helium atom scattering (HAS) studies of the H-covered Mo (110) and W(llO) surfaces reveal 
a twofold anomaly in the respective dispersion curves. In order to explain this unusual behavior 
we performed density functional theory calculations of the atomic and electronic structure, the 
vibrational properties, and the spectrum of electron-hole excitations of those surfaces. Our work 
provides evidence for hydrogen adsorption induced Fermi surface nesting. The respective nesting 
vectors are in excellent agreement with the HAS data and recent angle resolved photoemission 
experiments of the H-covered alloy system Moo.95Reo.05 (HO). Also, we investigated the electron- 
phonon coupling and discovered that the Rayleigh phonon frequency is lowered for those critical 
wave vectors. Moreover, the smaller indentation in the HAS spectra can be clearly identified as a 
Kohn anomaly. Based on our results for the susceptibility and the recently improved understanding 
of the He scattering mechanism we argue that the larger anomalous dip is due to a direct interaction 
of the He atoms with electron-hole excitations at the Fermi level. 

63.20.Kr, 68.35.Ja, 73.20.At, 73.20.Mf 



I. INTRODUCTION 



The interest in the (110) surfaces of Mo and W has 
been fostered in recent years by the discovery of deep and 
extremely sharp indentations in the energy loss speCitra of 
' He atom scattering at H/Mo (110) and H/W (110)BJ3. As 
I . depicted in Fig. ^ those anomalies are seen at an incom- 
. mensurate wavevector, Qg^pi along the [001] direction 

a (Pi?) and additionally at the commensurate wavevector 

^ ' Qcxp = 5" at the boundary of the surface Brillouin zone 
. ^ , (SBZ). At those points two simultaneous anomalies de- 
• velop out of the ordinary Rayleigh mode. One, wi, is 
^ I extremely deen._ajrujl is only seen by helium atom scat- 
tering (HAS)M'ffl. The other, uj2, is instead modest, 
and is observed by both HAS and high msdution elec- 
tron energy loss spectroscopy (HREELS)ErLl'El. Since the 

c,l 

[001] components of both critical wavevectors Q ' and 

c .2 

Q ' are approximately the same it was suggested that 
the anomaly runs parallel to PA^ through the SBZ as 
illustrated in Pig. ||. 

Various models have been brought forward in order 
to account for the unusual vibrational properties of H- 
covered Mo(llO) and W(llO): Por instance, a pro- 
nounced but less sharp softening of suiiace phonous 
is also seen on the (001) surfaces of WIIj and McO. 
There, the effect is cauaecLhji ^marked nesting proper- 
ties of the Permi surface^'EjO'Ej. However, angular re- 
solved photoemission (ARP) experiments of H/W (110) 
and H/Mo (110) gave no evidence for nesting vectors 
comparable to the HAS determined critical wave vec- 
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FIG. 1. HREELS spectra of the clean W(llO) surface (up- 
per panel) and its H-covered (1x1) phase (lower panel) from 
Ref. ^ Shown are the dispersion of the Rayleigh wave (trian- 
gles) and the longitudinally polarized surface phonons of the 
first (circles) and second (squares) layer. The dots reprint 
the results of Hulpke and Liidecke's HAS measurementsaa. 



TV P 




HAS anomalies 
[001] 



FIG. 2. Position of the HAS anomalies for the H/W(110) 
adsorbat system within the SBZ. The dotted line indicates 
the form of the anomaly apart from the symmetry directions 
VH and VS as suggested by Hulpke and Liideckeu. 
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FIG. 3. ARIiT-measured Fermi surface of the H-covered 
W(llO) surfacea. The shaded areas are the (110) projec- 
tion of the bulk Fermi surface. Surface states at the Fermi 
level are marked by dots and repeated in the other equivalent 
sections by crosses. The arrows represent the HAS critical 
wavevectors Q^'^^ and Q^^p- 

(see Fig. ^). Thus, for the phonon anoma- 
Ues of those two systems there seems to be no connec- 
tion between the electronic structure and the vibrational 
properties. Moreover, even if such a relation existed it 
would be unclear whether the phenomena has to be inter- 
preted as a Kohn anomaly or as the phason and ampli- 
tudon modes associated with the occurance of a charge- 
density wave. In fact, the experimental finding of two 
modes at S as well as some tliepretical arguments rule out 
the phason-amplitudon ideally. Furthermore, any model 
which links the jphonon anomalies to the motion of the 
hydrogen atomsQ has to be ruled out because the HAS 
spectra remain practically unchanged when deuterium is 
adsorbed instead of hydrogena^El. 

Still unexplained, but probably not directly correlated 
to the anomalies, is the peculiar behavior of hydrogen vi- 
braiions observed in HREELS experiments by Balden et 
al.uu. At a H-coverage of one mono-layer the adsorbate 
modes parallel to the surface form a continuum inter- 
preted a liquid-like phase of the adsorbate. 

An additional puzzle to the already confusing picture is 
added by the observation of a symmetry loss in the low- 
energy electron diffraction (LEED) pattern of W(llO) 
upon H-adsorption. The phenomena was believed being 
caused by a H-induced displaciancnt of the top layer W 
atoms along the [TlO] directiontj (see Figs.^ and c). For 
H/Mo(110) similar studies do notr^rovide any evidence 
for a top-layer-shift reconstructioncj. 

The goal of our work is to give a comprehensive expla- 
nation for the observed anomalous behavior and clarify 
the confusing picture drawn by the different experimen- 
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FIG. 4. (a) Adsorbate positions within the (1 x 1) unit 
cell of the W(llO) and Mo (110) surfaces: Shown are the 
long-bridge (L), short-bridge (S), hollow (H), and on-top (T) 
site, (b-c) Structure of-fche H-covered W(llO) surface as sug- 
gested by Chung et alx5: (b) low H-coverage 6 < 0.5 ML; (c) 
top-layer-shift reconstruction for a H-coverage > 0.5 ML. 
The white (shaded) circles represent the W atoms in the sur- 
face (subsurface) layer. The hydrogen positions are indicated 
by full dots, (d) Illustration of the structure parameters used 
in Table ^ The in-plane lattice constant is denoted by ao- 
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tal findings. We show that the anomalous behavior ob- 
served in the He and electron loss experiments is indeed 
governed by H-induced nesting features in the Fermi sur- 
face. 

The remainder of the paper is organized as follows. 
First, Section || gives a brief introduction into the ab ini- 
tio method used. The resu lts of o ur work are presented in 



Section [11: In Subsection III A we study the atomic ge- 



ometry of the clean and H-covered Mo(llO) and W(llO) 
surfaces. In particular, we check whether the hydrogen 
adsorption induces a top-layer-shift reconstruction. Fol- 



lowing in Subsection IIIB, we focus on the electronic 



structure of those surface systems. In this context, we are 
particularily interested in how the Fermi surface changes 
upon H-adsorption: Is there evidence for Fermi-surface 
nesting? If yes, how does it connect to the ARP measure- 
ments of the Fermi surfaces and the critical wavevectors 
detected by HAS and HREELS? Along the same line 



within the model of a Kohn anomaly, Subsection III C 
digs into the coupling between electronic surface states 
and surface phonons. There, we discuss the mechanism 



which causes the small anomaly cui . In Subsection [II D 
we examine the He scattering mechanism, study the spec- 
trum of electron-hole excitations, and explain the physi- 
cal nature of the huge anomalous branch uj2 ■ Finally, the 
major aspects of our work are summarized in Section 
The Appendix presents several detailed test calculations 
giving an insight into the properties and the accuracy of 
our method. 



II. METHOD 

All our calculations employ the density-functional the- 
ory (DFT) . They are performed within the framework of 
the local-density approximatiojiXLDA) for the exchange- 
correlation energy functionaSO. For the self-consistent 
solution of the Kohn-Sham (KS) equations we employ 
a full-potential . lijae arized augmented plane- wave (FP- 
LAPW) codeE3'ESE3 which_,we enhanced by the direct 



calculation of atomic forcesE3't3. Combined with damped 
Newton dynamicsES this enables an efhcient determina- 
tion of fully relaxed atomic geometries. Also, it allows 
the fast evaluation of phonon frequencies within a frozen- 
phonon approach. 



A. LAPW Formalism 

Thece|--je^sts_,a--j^^ different LAPW pro- 

gramJS^lHaayHaNeElElH. We therefore briefly 
specify the main features of our code. This is also in 
order to clarify the meaning of the LAPW parameters 
which determine the numerical accuracy of the calcula- 
tions. 

In the augmented plane-wave method space is di- 
vided into the interstitial region (IR) and non-overlap- 
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pingj— muffiii-tin (MT) spheres centered at the atomic 
sitesE3o'EJ. This division accounts for the atomic-hke 
character of the wavefunctions, the potential, and the 
electron density close to the nuclei and that the behav- 
ior of these quantities is smoother in between the atoms. 
The basis functions 0k (r) which we employ in our code 
for the expansion of the electron wavefunctions of the KS 
equation 

V'M(r)= C'.(K)0K(r) (1) 

|K|<ic:~f 



are defined as 



f7-i/2cxp(iKr), reIR 



0K(r)=<j Y,[aUK)uiiei,rj) + bUT^)ui{ei,rj)]YUri), rj<.sj . (2) 

Irn 

Here, K = k + G denotes the sum of a reciprocal lat- 
tice vector G and a vector k within the first Brillouin 
zone. The wavefunction cutoff K""^ limits the number of 
those K vectors and thus the size of the basis set. The 
other symbols in eq. have the following meaning: ft 
is the unit cell volume, s/ specifies the MT radius, and 
r/ = r — R/ represents a vector within the MT sphere of 
the I-th atom. Note that Yim{f) is a complex spherical 
harmonic with y;_„j(f) = (~^)"^Yi^{-f). The radial func- 
tions ui{ei,r) and ui{ei,r) are solutions of the equations 

H^P''ui{ei,r)Yi^{f) = q r)rj„(f ) (3) 
H'P''ui{ei,r)YUf) = [qu;(e,, r) + r)]l,™(f) (4) 

which are regular at the origin. Here, the KS operator 
^sph contains only the spherical average of the effective 
potential within the respective MT. The expansion ener- 
gies ei are chosen somewhere within the respective energy 
bands with Z-character, whereas the coefficients aj^^lK.) 
and 6^,„(K) are fixed by requiring that value and slope 
of the basis functions are continuous at the surface of the 
I-th MT sphere. Obviously, this can only be fulfilled for 
I < r^, but as is typically high (e.g. > 8) this 
does not represent a problem. 

The representation of the potentials and densities re- 
sembles the one employed for the wavefunctions, i.e., 

' VGexp(iGr), r e IR 

G|<Gp°' 

/ , Vim.i{ri)Yijn{ri), ri < Si . 

Im 



V{r) = < 



Thus, no shape approximation is introduced. The qual- 
ity of this full-potential description is controlled by the 
cutoff parameter Gp°* for the reciprocal lattice vectors G 
and the maximum angular momentum of the (Z, m)- 
representation inside the MTs. 
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B. Slab Systems and Parametrization 



The substrate surfaces are modeled by five-, seven-, 
and nine-layer slabs repeated periodically and separated 
by a vacuum region whose thickness is equivalent to four 
substrate layers. The MT radii for the W and Mo atoms 
are chosen to be smo = svv = 1.27 A. Note that the 
experimental interatomic distances of bulk Mo and W 
are 2.73 Aand 2.74 A, respectively. For hydrogen the MT 
radius is set to sh = 0.48 A. 

In the case of the W surfaces the valence and semi-core 
electrons are treated scalar-relativistically while the core 
electrons are handled fully relativistically. For Mo all 
calculations are done non-relativistically. In both met- 
als we introduce a second energy window for the treat- 
ment of the semi-core electrons (4s and 4p for Mo, 5s 
and 5p for W). The in-plane lattice constants (compare 
Fig. |) a^4f° = 3.14 A and a*^^™ = 3.13 A are calculated 
to be without including zero-point vibrations. They are 
in good agreement with the respective measured bulk 
lattice parameters at room t emp erature and other theo- 
retical results (see Appendix A 1 ). 

For the potential the (Z, m) representation within 
each MT sphere is taken up to = 4 while the 

kinetic-energy cutoff for the interstitial region is set to 
^(•(^wf-)2 ^ 100 Ry. Generally, we choose the plane- 
wave cutoff for the wavefunctions to be ^(K'"^)'^ = 
12 Ry; only for the frozen-phonon calculations we use 
a shghtly smaller cutoff of ^{K^^)^ = lORy which is 
indeed sufficient. The {I, m) representation of the wave- 
functions within the MTs is taken up to V"^ = 8. Also, 
we employ Fermi smearing with a broadening of k^T^i = 
68 meV in order to stabilize the self-consistency and the 
kj| -summation. All energies given below are however ex- 
trapolations to zero temperature Tei^^- We performed 
systematic tests comparias-the LDA and the generalized 
gradient approximationcfla exchange-correlation func- 
tional. For the quantities reported in this paper (to- 
tal energy differences of different adsorption sites, bond 
lengths, etc.) both treatments give practically the same 
results. All k sets used were determined by using the spe- 
cial k point scheme of Monkhorst and Fackea. Also, we 
checked the k|| -point convergence: In the case of atomic 
and electronic structure calculations a two-dimensional 
uniform mesh of 64 ky-points within the (1 x 1) SBZ 
gives stable results. For the evaluation of frozen-phonon 
energies a set of 56 k|| -points is employed within the SBZ 
of th e enla rged (1x2) and (2x1) surface cells (see subsec- 
tion A 3 b ) . In Table | we summarize the main features 
of our LAPW parameter setting. For more details the 
reader is referred to the discussion in the Appendix. 
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TABLE I. LAPW parameter setting used in Section |D| for 
the bulk studies, the calculation of the atomic and electronic 
surface strucutures, and the evaluation of frozen-phonons. 



LAPW 
parameter 


bulk 


surface 
structure phonons 


pot 


4 


4 


3 


^^(GP-)2 [Ry] 
&iK-r [Ry] 


64 
8 


100 
8 


100 
8 


12 


12 


10 


Fermi smearing [meV] 


68 


68 


68 


number of k-points 








for valence electrons 


3375 


64 


56 


for semicore electrons 


216 


9 


4 


in unit cell 


bulk bcc 


(ixl) 


(1x2) 








(2x1) 



III. RESULTS 



A. Atomic Structure 



Using the directly calculated forces in combination 
with a damped Newton dynamics of the nuclei the atomic 
structure of stable and metastable geometries is obtained 
quite automatically. It is important that the starting con- 
figuration for an optimization run is of low (or no) sym- 
metry. Otherwise the system may not relax into a low- 
symmetry structure of possibly lower energy. Also, inde- 
pendent of the optimization method used, there is always 
a chance that the minimization of the total energy leads 
into a local minimum instead of the global one. In order 
to reduce and hopefully abolish this risk one has to con- 
duct several optimizations with strongly varying starting 
configurations. We consider a system to be in a stable (or 
meta-stable) geometry if all force components are smaller 
than 20 me V/A. The error in the structure parameters of 
the relaxed system is then ±0.02 A « ±1.0%do for the W 
systems a nd ± 0.01 A w ±0.5%(io for the Mo surfaces (see 
Appendix A 2). 

The relaxation parameters calculated for the clean and 
H-covered (110) surfaces of Mo and W are presented in 
Table |l| (Fig. |jd serves to clarify the meaning of these 
parameters). The converged nine-layer-slab results are 
remarkably similar for both transition metals. Moreover, 
they show excellent quantitative agreement with the re- 
sults of a recent LEED analysis also presented in Ta- 
ble |1. For a detailed comparision between experiment 
and theory we refer to Ref. ^ As the energetically most 
favorable hydrogen adsorption site on both W (110) and 
Mo (110) we identify a quasi-threefold position (indicated 
as "H" in Fig. ^). The calculated adsorption energies 
for other possible positions shown in Fig. ^ are several 
100 meV less favorabldfj. 

Our investigations also throw light on the suggested 
model of a H-induced structural change: For both mate- 
rials the calculated shift yi is only of the order of 0.01 A 
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TABLE II. Calculated relaxation parameters for the clean and H-covered (110) surfaces of Mo and W. The height of the 
hydrogen above the surface and its [TlO] offset from the [001] bridge position are denoted by dn and j/h, respectively. The shift 
of the surface layer with respect to the substrate is j/i. The parameters Adij describe the percentage change of the interlayer 
distance between the i-th and the j'-th substrate layers with respect to the bulk interlaypfiiseacing do- For each system the 
results of a five- , seven- , and nine-layer calculation as well as of a recent LEED analysifc3c3 are presented (labelled as "5", 
"7", "9", and "L"). 



system 






dn 


Vi- 


Adi2 


Ad23 


Ad34 










rXi 


[/odo\ 


[/odo\ 


roy J 1 
[/odo\ 


Mo(llO) 


5 








—5.9 


—0.8 






7 








-4.5 


±0.5 


0.0 




9 








-5.0 


±0.7 


—0.3 




L 








-4.0 ±0.6 


±0.2 ±0.8 


0.0 ± 1.1 


H/Mo(110) 


5 


0.63 


1.08 


0.05 


-2.7 


-0.4 






7 


0.60 


1.07 


0.03 


-2.1 


±0.1 


-0.1 




9 


0.62 


1.09 


0.04 


-2.5 


±0.3 


-0.2 




L 


0.55 ±0.4 


1.3 ±0.3 


0.0 ±0.1 


-2.0 ±0.4 


0.0 ±0.5 


0.0 ±0.8 


W(llO) 


5 








-4.1 


-0.2 






7 








-3.3 


-0.1 


-0.4 




9 








-3.6 


±0.2 


-0.3 




L 








-3.1 ±0.6 


0.0 ±0.9 


0.0 ± 1.0 


H/W(110) 


5 


0.68 


1.12 


0.01 


-1.4 


0.0 






7 


0.67 


1.11 


0.05 


-1.3 


±0.3 


±0.3 




9 


0.67 


1.09 


0.02 


-1.4 


-0.3 


-0.1 




L 


0.56 ±0.4 


1.2 ±0.25 


0.0 ±0.1 


-1.7 ±0.5 


0.0 ±0.6 


0.0±0.9 




FIG. 5. Change of total energy _E*° versus top-layer-shift 
yi calculated for a H/W(110) seven-layer slab system. For 
each data point the whole surface is relaxed keeping only the 
substrate [110]-coordinates yi and 2/2 = 3/3 = A fixed. Also 
shown are the first three oscillator eigenstates calculated from 
a harmonic expansion of the total energy _E'°'(t/i). 
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and thus there is no evidence for a pronounced top-layer- 
shift reconstruction. Moreover, this subtle change in the 
surface geometry is unlikely to be resolved experimen- 
tally due to zero-point vibrations. This aspect becomes 
evident by evaluating the total energy with respect to 
a rigid top-layer-shift. In Fig. [s] we present the results 
of such a calculation for H/W(110) and depict the first 
three vibrational eigenstates obtained from a harmonic 
expansion of the total energy i?*°*(yi); for H/Mo(110) 
the energetics is similar. Since we have k^T « 25meV 
at room temperature thermal fluctuations of j/i are of the 
order of 0.1 A. This is considerably larger than the the- 
oretically predicted ground state value of yi. In Ref. ^ 
Arnold et al. offer an explanation of how the previous 
LEED work was probably mislead by assuming that scat- 
tering from the H layer is negligible. 

B. Electronic Structure 

In Section | we already pointed out that one needs to 
find pronounced nesting features in the Fermi surface of 
the H-covered surfaces in order to make the model of a 
Kohn anomaly work. After the determination of the re- 
laxed surface configurations we therefore turn our atten- 
tion towards the electronic structure of those systems. In 
particular, we focus on surface states close to the Fermi 
level. 

One should keep in mind that DFT is not expected to 
give the exact Fermi surface, i.e., the self-energy operator 
may have (a strong or weak) k depencence, different from 
that of the exchange-correlation potential. Moreover, we 
encounter additional problems because it is rather diffi- 
cult in slab systems to distinguish between surface reso- 
nances and pure surface states. Due to interactions be- 
tween the two faces on either side of the slab resonances 
can be shifted into a bulk band gap. There, they are eas- 
ily mixed up with real surface states. By contrast, in a 
semi-infinite substrate only true surface states are local- 
ized in the bulk band gap and decay exponentially into 
the bulk. Therefore, our slab method lacks a clear indi- 
cator for the characterization of electronic states. How- 
ever, we found that the localization w^'^ within the MT 
spheres of the top two surface layers of a seven-layer slab 
is a useful and suitable measure (60% < w^'^: surface- 
state like, 30% < w^'^ < 60%: surface -resonance like, 
^MT ^ 30%: bulk like). Nevertheless, this approach is 
anything but exact. Also, the calculated k|| space loca- 
tion is only accurate if the respective band is strongly 
localized at the surface. 

Experimentalists face similarly serious problems: In 
ARP the Fermi surface is determined by extrapolation 
peaks found below the Fermi level. However, the value 
of the Fermi energy itself is only known within an uncer- 
tainty of about ±100 meV. This uncertainty can amount 
to a significant error in the extrapolated k|| space position 
of states, especially if the respective bands are flatE3. 
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theory experiment 




FIG. 6. Theoretical and ARP Fermi surfaces of the clean 
(upper part) und the H-covered (lower part) W(llO) surface. 
The solid (dotted) lines denote surface resonances or surface 
states which are localized by more than 60 % (30 %) in the 
MTs of the two top W layers. Shaded areas represent the 
(110) projected theoretical W bulk Fermi surface. The arrows 

c 1 c 2 

Qcxp ^^'^ Q ''■'"^ tlis critical wave vectors found in HAS. The 
ARP data stem from Refs. ^, |l^. 
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theory experiment 




FIG. 7. Theoretical Fermi surfaces of the (a) clean and (c) 
H-covered Mo (110) surface. Also presented are data points 
(+) which stem from an ARP study le^ the (b) clean and (d) 
H-covered Moo.95Reo.05 (HO) surfaceEI. The presentation is 
equivalent to Fig. 0. 
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FIG. 8. Hydrogen induced shift of the {d^ 



, dxz) band 



at the Mo(llO) surface. The upper part of the figure rep- 
resents a one-dimensional cut through the two-dimensional 
surface bandstructure. Shaded areas illustrate the (110) pro- 
jected bulk band structure. The dotted (solid) lines depict the 
position of the {d^^2_^2,dxz) band for the clean (H-covered) 
surface. The arrows indicate how the bandstructure trans- 
lates into the two-dimensional Fermi surface shown in the 
lower part. In Fig. |^ the charge density of the two states 
marked by the arrows is presented. 



In Figs. 1^ and we report the calculated Fermi sur- 
faces and Gopjpare them to expf;rimental ones obtained 
by Kevan'sEl'EII and Plummer'sEEl groups. Let us first fo- 
cus on the theoretical results which are — as in the case 
of the surface geometries — again very similar for Mo 
and W. 

For both systems the H-adsorption induces the shift 
of a bandAwith (d^z^-r^jd^z) character to lower binding 
energiesEjO. This effect which is illustrated in Fig. g 
moves the Fermi line associated with this band into 
the band gap of the surface projected band structure. 
Subsequently, the respective states become true surface 
states. It is important to note that the band shift is 
due not to a hybridization between hydrogen orbitals and 
{d3z2_j.2,dxz) bands but to a hydrogen induced modifi- 
cation of the surface potential. The bonding states of 
the hydrogen-substrate interaction are about 5 eV below 
the Fermi energy and the anti-bonding states are 4eV 
above. Therefore, they are not involved in this process 
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Mo(llO) 



H/Mo(110) 



[10-VBohr'] 




FIG. 9. Charge density plot along the (010) plane for 
Mo(llO) and H/Mo(110): Shown are the two states marked in 
Fig. ^ by arrows. Only the plane-wave part of the wavefunc- 
tions is shown. The positions of the Mo atoms are marked by 
dots. The square represents the hydrogen atom in the hollow 
adsorption site. 



which takes place at the Fermi level. 

The shifted {d^z^_r2,dxz) band is characterized by a 
high density of states at the Fermi level. For the clean 
Mo(llO) surface one finds that the {d3^2_^2^dxz) band 
has a MT locahzation of w^^ « 30%. Due to the H- 
adsorption this value increases to more than 60%. This 
modification is also visible in the charge density plots 
in Fig. ^. More important, the new Fermi contour re- 
veals pronounced nesting features. In Figs, ^c and Qc one 
sees the dramatic changes due to the H-induced shift. 
Segments of the two-dimensional Fermi contour of the 
{dsz^-r^ , dxz) band are now running parallel to TN and 
normal to TS. There are two nesting vectors Q"^'^ and 



TABLE III. Theoretical Fermi surface nesting vectors com- 
pared to critip^iyavevectors obtained by HAS and HREELS 
experimentaJt 



direction 


system 










theory 


experiment 


FH 


H/Mo(110) 


0.86 


0.90 




H/W(110) 


0.96 


0.95 


FS 


H/Mo(110) 


1.23 


1.22 




H/W(110) 


1.22 


1.22 
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c 2 . 

Q ' which connect those segments in different parts of 
the SBZ. As can be seen in Table ID they agree very well 
with the HAS measured critical wavevectors. 

With respect to those nesting vectors our results con- 
tradk±_Lh£, photoemission studies by Kevan and cowork- 
eraj^EElEZIH which were depicted in Figs. ||b,d (for W) and 
Fig. 2 in Ref. ^ (for Mo). Apart from that Kevan's and 
our data compare rather well. Therefore, it is difficult to 
give a plausible reason for the discrepancies. However, in 
view of the fact that for the clean surface the theoretical 
and experimental Fermi surfaces are in very good agree- 
ment and that our calculated Fermi nesting vectors agree 
very well with the HAS and HREELS anomalies we dared 
to suggest that those differences may be due to the al- 
readjy discussed problems within the experimental analy- 
sisC3. We also note that recent ARP studied^ of a related 
system which we present in Fig. ^ seem to support our 
conclusion. This experimental work deals with the (110) 
surface of the alloy Moo.95Reo.05 but the surface physics 
of Moo.95Reo.05 (110) and Mo (110) should be practically 
the same because in both cases the top layer consists 
only of Mo atoms. This assumption is also backed by 
test calculations where we simulated a Moo.95Tco.05 al- 
loy be the virtual crystal approximation and found only 
minor quantitative changes. 

From Figs. 0c and 0d it becomes clear that in particu- 
lar for the important (d^z^-r^^dxz) surface band, which 
was not seen by Kevan's group, experiment and DFT 
now agree very well. There are however still differences: 
Theory predicts bands centered at S which are not seen 
by ARP whereas the band circle at F in Fig. ^ is only 
observed experimentally. Also, in the calculations we find 
an elliptical band centered at the F point whereas in ARP 
the same band has the shape of a rectangle. Those dis- 
crepancies are probably due to matrix elements of the 
photoemission process and to the inherent inaccuracies in 
theory and experiment which we mentioned above. Nev- 
ertheless, the theoretical results are encouraging because 
they provide strong evidence for a link between the HAS 
anomalies and the Fermi surfaces. 



C. Vibrational Properties 



The study of the electronic structure revealed pro- 
nounced nesting features for the H-covered surfaces. At 
this point, the following question arises: How does the 
surface react to this apparent electronic instability? At 
the k vectors of the Fermi surface nesting the coupling 
between electrons and phonons is expected to become sig- 
nificant which implies a possible breakdown of the Born- 
Oppenheimer approximation. This leads to a softening 
of the related phonons. If the electron-phonon coupling 
is strong and the energetic cost of a surface distortion 
is small the nesting could even trigger a reconstruction 
combined with the build-up of a chaxge-densitK-wave as 
in the case of the (001) surfaces of WEI and MoB. Then, 
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at the reconstructed surface, the Born-Oppenheimer ap- 
proximation is vaUd again. It is clear that one needs to 
perform frozen-phonon calculations in order to determine 
the actual strength of the electtpn-phonon coupling and 
the resulting phonon softcningtS. 

The experimental vibrational spectra of the (110) sur- 
faces of Mo and W show two distinct reactions to the 
adsorption of hydrogen. Along TH and TS the H- 
adsorption induces a softening of the RayJfiigh and (to a 
smaller amount) of the longitudinal waveEj while a stiff- 
ening of these modes is observed along TN. Our goal is 
to investigate both effects theoretically. 

We use an enlarged surface unit cell together with a 
five layer slab. The plane-wave cutoff for the wavefunc- 
tions is reduced to iL- (j^' ^f)2 — 10 Ry (see test cacula- 
tions in the Appendix |A 5| ) . The geometries are defined 
by the five layer relaxation parameters presented in Ta- 
ble Some relaxation parameters, i.e., Ady, change 
considerably when we perform a calculation with a slab 
of five instead of seven or nine metal layers. We note, 
however, that the values of the critical wavevectors and 
the position of the hydrogen with respect to the substrate 
surface are practically insensitive to the thickness of the 
slab. This indicates that the physical properties we are 
interested in, e.g. the nesting features, are well localized 
surface phenomena and that the results of our five layer 
slab studies can be trusted. As shown in Appendix A3 
for the iS* point phonon the calculated frequencies are rel- 
atively insensitive to the the SBZ k|| point sampling, and 
we conclude that a uniform mesh of 56 points is sufficient 
for the present study. 

Being zone boundary phonons in the Rayleigh waves 
at S and N the displacements of the surface atoms are 
along the direction normal to the surface. Moreover, the 
vibrations are strongly localized in the surface region. 
We assume that even with hydrogen adsorbed the cou- 
pling to modes parallel to the surface is small. Therefore, 
we can confine our study to the vibrational components 
normal to the surface. Besides the first substrate layer 
we also include the vibrations of the second layer in our 
calculations. Due to the small mass the hydrogen atom 
follows the substrate vibrations adiabatically. In the cal- 
culations this adsorbate relaxation lowers the distortion 
energy by several meV and thus reduces the phonon en- 
ergy significantly. Therefore, all phonon frequencies are 
calculated for a fully relaxed hydrogen adsorbate. 

For the study of the vibrational properties of the sys- 
tems we displace the substrate atoms in the surface (ac- 
cording to Fig. |lO|), relax the hydrogen atom, and calcu- 
late the resulting atomic forces. The same is done for the 
atoms in the subsurface layer. From a third-order fit to 
the calculated forces for five different displacement steps 
between 0% and 5% of the lattice constant we obtain the 
matrix elements of the dynamical matrix and then via 
diagonalization the respective phonon energies. The cal- 
culated frequencies are collected in Table IV. At the N 
point our results reproduce the experimentally observed 
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FIG. 10. Distortion pattern of the atoms of the top metal 
layer for the Rayleigh phonons at the symmetry points S and 
iV. 

TABLE IV. Comparison of calculated frozen-pho«pji en- 
ergies and experimental values obtained by HASolI and 
HREELSa. The theoretical phonon energies are obtained us- 
ing a five-layer slab. Their numerical accuracy is about ±5%, 
i.e., ±1 meV. 



phonon 


system 




[meV] 






theory 


experiment 


N 


W(llO) 


15.4 


14.5 




H/W(110) 


17.6 


17.0 


S 


Mo (110) 


22.7 


~21 




H/Mo(110) 


17.2 


<16 




W(llO) 


18.3 


16.1 




H/W(110) 


12.0 


11.0 



increase of the Rayleigh-virave frequency as hydrogen is 
adsorbed. 

At the S point we find that the strong coupling to 
electronic states at the Fermi level leads to a lowering of 
the phonon energy in good agreement with the experi- 
mental results. In Figs. pl]a-c we schematically illustrate 
the mechanism which is responsible for this effect. It 
is, in fact, a text-book epapiple of a Kohn anomaly due 
to Fermi-surface nestingEJ'E3: (a) Within the unpertur- 
bated system the {d^z'^-r^, d^z) band cuts the Fermi level 
exactly midway between F and S. (b-c) The nuclear dis- 
tortion associated with the S-point phonon modifies the 
surface potential and hence removes the degeneracy at 
the backfolded zone boundary point S'. The occupied 
states are shifted to lower energies (see Fig. ^2|). This 
amounts to a negative contribution of the electronic band 
structure energy to the total energy and thus a lowering 
of the phonon energy. 

A frozen-phonon study for the second nesting vector 
along TH is not performed because the respective Q'^' 
is highly non-commensurate. Thus, such a calculation 
would be very expensive. However, since the character 
of the {d^z^_j.2, d^z) band does not change when shift- 
ing from the TS to the TH nesting we exped— similar 
results; this was recently confirmed by BungaroEj within 
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FIG. 11. Schematic representation of the mechanism which 
causes the Kohn anomaly of the H/Mo(110) and H/W(110). 
Shown is the bandstructure along SFS parallel to the nesting 
vector Q'^' in Figs. ^ and |^. The form of the {d^z^-r^ i dxz) 
band is indicated by dashed lines. See text for a detailed 
description. 
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FIG. 12. H/W(110): Band gap A induced by a 5 point 
frozen-phonon distortion. Shown are the KS eigenvalues of 
the {d^z^-r^ J dxz) band at the S' point versus the amplitude 
of the atomic distortion tiz of the first substrate layer. 
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the framework of a perturbation theory. 

The calculation of the electron-phonon interaction at 
the Mo (110) and W (110) surfaces and its change due to 
hydrogen adsorption pinpoints the phonon character of 
the small anomaly L02 and identifies the interplay between 
the electronic structure and the vibrational spectra of 
the transition metal surfaces. Thus, our results clearly 
support the interpretation that the small dip observed 
by both HAS and HREELS is due to a Kohn anomaly. 
Furthermore, we find that the electron-phonon coupling 
is not strong enough to induce a stable reconstruction. 
Thus, the system remains in a somewhat limbo-like state. 

D. Electron- Hole Excitations 

For the deep and narrow anomaly lOx the interpretation 
appears to be less straightforward. It is particularly puz- 
zling that it is only seen in the HAS spectra and not in 
HREELS. This difference, together with the above noted 
"limbo state" of the surface yield the clue to our present 
interpretation. At first, it is necessary to understand the 
nature of rare-gas atom scattering. 

In a recent study we found that those scattering pro- 
cesses are significantly more coroalicated (and more in- 
teresting) than hitherto assumedE2: The reflection of the 
He atom happens in front of the surface at a distance 
of 2-3 A. This is illustrated in Fig. |l^. More important, 
it is not the total electron density of the substrate sur- 
face which determines the interaction but the electronic 
wavefunctions close to the Fermi level. In the case of 
the H/W(110) and H/Mo(110) systems it is thus plau- 
sible to assume that the He atom couples directly to 
the ((i3z2_r2 , dxz) surface states mentioned above and ex- 
cites electron- hole pairs. By contrast, the electrons in 
HREELS scatter at the atomic cores and interact only 
weakly with the electron density at the surface. 

In order to analyze the spectrum of those excitations 
as seen by HAS in some more detail we evaluate the local 
density of electron-hole excitations 

P(q||,?iw)= / dk|| u^kii+qiiU'kii (/kn+qn -/k,,) 

JSBZ 

X(5(ek||+qj| - ek|| - txijj) (6) 

which is a measure of the probability that a He atom 
looses the energy Tilo and the momentum q|| when inter- 
acting with the surface elctron density. In our approach 
we use all eigenvalues eky and occupation numbers /ky 
obtained via a nine-layer slab calculations. In order to 
refer to HAS we take into account the localization Wky 
of the respective state at a distance of 2.5 A in front of 
the substrate surface. 

The results obtained for Mo (110) and H/Mo (110) are 
presented in Figs. For the clean surface we find that 
the intensity of the electron-hole excitations decreases 
continuously as we move away from the F point towards 
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He atom scattering 




Electron scattering 




FIG. 13. Schematic illustration of the scattering of He 
atoms (HAS, upper part) and electrons (HREELS, lower part) 
at a metal surface (from Ref. M). 
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FIG. 14. Contour plot of the local probability func- 
tion P(q||,Rcj) of Mo (110) and H/Mo(110) calculated for 
Jko = 0.5 mRy « 6.8 meV. The dashed line represents a value 
of 3 arbitrary units while each full line denotes an additional 
increase by 1 arbitrary units (with increasing line width). The 
positions of the HAS measured anomalies within the SBZ are 
indicated by dots. 
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W(llO) 




the zone boundaries. This relatively smooth behavior 
of the susceptibility is modified considerably when hy- 
drogen is adsorbed on the clean surface: Pronounced 
peaks appear, in particular one which stretches from 
the TH line through the SBZ to the symmetry point 
S. This is in excellent agreement with the HAS results 
(see also Fig. ||). Also, there is a peak located at the 
upper half of the TN symmetry line. It cannot be as- 
sociated with an anomaly in the HAS spectrum. This 
result may be due to the fact that selection rules or ma- 
trix elements for the interaction between the scattering 
He atoms and the electron-hole excitations are not con- 
sidered in P{'k,huj). In Fig. |l^ the respective data for 
the clean and H-covered W(llO) surfaces are displayed. 
In general, our findings are in agreement with the exper- 
imental results of Liidecke and Hulpke that the anomaly 
stretches through the SBZ. However, we do not support 
the idea of a fully one-dimensional nesting mechanism 
with a constant [001]-component of the critical wave vec- 
tor as depicted in Fig. ^. 

In view of the simplicity of our approach, the calcu- 
lated spectrum of electron-hole excitations as predicted 
by P(q, hu) seems to be well-connected to the HAS mea- 
surements. Moreover, it is surprising how detailed the 
Fermi-surface nesting (compare to Figs. ^ and 0c) is 
still present at a (rather large) distance from the surface. 
In fact, if the localization were not included into 
the calculations or if it were taken closer to the surface 
the H-induced peaks would be hidden in the background 
noise of other electron-hole excitations. The clear image 
in Figs. |l^ and |l^ only appears if we take the position of 
the classical turning point into account. 

To our knowledge these are the first theoretical studies 
of the rare-gas atom scattering spectrum using a first- 
principles electronic bandstructure. Our findings support 
the earlier suggestionCJ that the giant indentations, seen 
only in HAS, are due to the excitation of electron-hole 
pairs. In HREELS such excitations are less likely, and 
thus the strong anomaly remains invisible. 

IV. SUMMARY AND OUTLOOK 

In conclusion, the major results of our work are as 
follows: We demonstrate that for both W(llO) and 
Mo (110) the (1 X 1) surface geometry is stable upon 
hydrogen adsorption. There is no evidence for a pro- 
nounced H-induced top-layer-shift reconstatction, a re- 
sult confirmed by a recent LEED analysisc3'C3. Also, we 
give a consistent explanation for the H-induced anoma- 
lies in the HAS spectra of W (110) and Mo (110). The H- 
adsorption induces surface states of {d^y2_j.2^dxz) char- 
acter that show pronounced Fermi-surface nesting. The 
modestly deep anomaly is idpatiJled as a Kohn anomaly 
due to those nesting featureaij'El'Ea. For an understand- 
ing of the huge dip we stress that the scattering of rare- 
gas atoms is crucially influenced by interactions with 
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substrate|-S|Urface wavefunctions at the classical turn- 
ing pointEa. Assuming that the He atom couples effi- 
ciently to the H-induced surface states on H/W (110) and 
H/Mo(110) we therefore conclude that the deep HAS 
anomaly is predominantly caused by a direct excitation 
of electron-hole pairs during the scattering process. 

We hope that our calculations and interpretations 
stimulate additional theoretical and experimental work: 
For instance, it is highly important to better understand 
the details of rare-gas atom scattering processes at "real" 
surfaces. Also, the liquid-like behavior of the H-adatoms 
observed in the HREELS measurcmentsa remains an un- 
resolved puzzle. One would also like to know whether 
the HAS and HREELS spectra of the Moo.95Reo.05 (HO) 
alloy system reveal H-induccd anomalies as expected. Fi- 
nally, we call for a new ARP study of the clean and H- 
covered Mo (110) and W (110) surfaces in order to exper- 
imentally identify the Fermi surface. We also note that 
scattering experiments with atoms or molecules like Ne 
or H2 might provide additional insight into the interest- 
ing behavior of these surfaces. 
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APPENDIX: TEST CALCULATIONS 
1. Bulk 

In Table we list the calculated equilibrium lattice 
constants qq and the bulk moduH Bq = ~V d'^ E*-°^ / d^V , 
and compare them to published results from experiment 
and theory. Our values were obtained using the LAPW 
parameter setting presented in Table |[ 

2. Surfaces 

a. Test System 

The focus of this paper lies on the Mo(llO) and W(llO) 
surfaces. In order to test the numerical accuracy of our 
calculations for such systems it is important to go beyond 
typical tests for the crystal bulk. For instance, we need 
to know the error bars for atomic geometries obtained 
via forces as compared to a pure total energy calcula- 
tion. Five- and more-layer systems are not suitable and 
in fact for most cases also not necessary for such a de- 
tailed study. Therefore, questions refering to the size of 
the vacuum region, to the accuracy of the directly calcu- 
lated LAPW forces, and to the appropriate wavefunction 



24 



TABLE V. Theoretically and experimentally obtained crystal parameters for Mo and W. The 
percentage deviations from experiment are given in parentheses. The results of this work are 
calculated non-relativistically (N) and scalar-rekiitivistically (R). Also, tha-different XC potentials 
used in the-theoretical studies are given: WignerEj (W37), Barth and HedinESI (BH72), and Ceperley 
and Aldeitil (CA80) parameterized by Perdew and ZungeiEa. The abreviation "NL-PP" stands for 
"non-local pseudo-potential" method. All theoretical numbers ignore the influence of zero-point 
vibrations. 



metal 


method 


XC 


ao [A] 




Bo [MBar] 


Mo 


experiment^ 




3.148 




2.608 




NL-PP*" 


BH72 


3.152 


(+0.1) 


3.05 (15.5) 




LAPW^ 


BH72 


3.131 


(-0.5) 


2.91 (11.6) 




this work (N) 


CA80 


3.126 


(-0.7) 


2.88 (10.4) 




this work (R) 


CA80 


3.115 


(-1.1) 


2.89 (10.8) 


W 


experiment'' 




3.163 




3.23 




NL-PP*" 


BH72 


3.173 


(+0.3) 


3.45 (6.8) 




LAPW" 


W34 


3.149 


(-0.4) 


3.46 (7.1) 




LAPW= 


BH72 


3.162 


(±0.0) 


3.40 (5.3) 




this work (N) 


CA80 


3.194 


(+1.0) 


2.92 (-^9.6) 




this work (R) 


CA80 


3.137 


(-0.8) 


3.37 (4.3) 


^from Ref. 



'"from Ref. H 
''from Ref. PSf 
''from Refsr|63, || 
'from Ref. [34| 



cutoff will be answered by studying a simple Mo(llO) 
double-layer slab system. Starting from the LAPW pa- 
rameter set listed in Table | we vary one parameter and 
study its influence on the evaluated total energy and the 
atomic force. These quantities are calculated with re- 
spect to the inter-layer distance z. Again, we note that 
we are interested not in the physics of this double-layer 
system but in the accuracy and the convergence of our 
method. Therefore, the calculations are done using just 
four k|| points within the SBZ. 

Employing a Jaast-square fit which has the form of a 
Morse potentiaE3 

E''°\z) = d\\-c-^^''-'^~'^'^^'^ (Al) 

we calculate the equilibrium relaxation parameter d. 
Here, c?o represents (110) layer distance in the bcc crystal. 
Additionally, the vibrational properties are of consider- 
able interest. From the Morse parameters D and /3 one 
extracts the oscillator frequency by using the expression 

where the atomic masses M/ are those of the Mo and 
W nuclei. A corresponding fit is also employed for the 
directly calculated forces normal to the (110) surface: 

dz 

= — 2Z)/3 Q-'^Mz-d-do) _ ^-f3(z-d-do) ^ {^i) 

This enables a quantitative study of the accuracy of the 
evaluated LAPW forces with respect to the total energy 
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TABLE VI. Mo(llO) double-layer: Morse parameter d and 
the energy Tuu versus the number of vacuum layers. The arrow 
indicates the parameter choice which yields converged results. 



number d [%cio] Tiuj [meV] 



of layers 


E 


F 


E 


F 


2 


-1.7 


-2.1 


36.5 


35.3 


3 


-1.7 


-2.1 


36.7 


35.2 


=^ 4 


-2.0 


-2.4 


37.3 


35.6 


5 


-2.1 


-2.4 


37.4 


35.7 


6 


-2.0 


-2.3 


36.9 


35.4 



TABLE VII. Mo(llO) double-layer: Variation of the wave- 
function cutoff parameter ^(^^"')^. The notation is equiv- 
alent to Table fvil 



[Ry] 


d [%do] 




Tiuj 


[moV] 


E 


F 


E 


F 


9 


-0.5 


-1.0 


35.0 


33.4 


10 


-1.5 


-2.1 


36.0 


34.6 


11 


-1.6 


-2.2 


36.1 


35.2 


=^ 12 


-2.0 


-2.4 


37.3 


35.6 


13 


-1.9 


-2.2 


37.1 


35.5 


14 


-1.9 


-2.1 


36.5 


35.4 



E^°\z). Within the Tables 
are marked by "F" (force) and "E" (energy). 

We consider a calculation to be converged if the re- 
laxation parameter d and the frequency oj differ by less 
the Ad = 0.25%^^=°"^ « 0.005 A and ^iAw = ;i5%w™"^ « 
2meV from the fully converged results d'^™^ and . 
In the following discussion we focus on these two quan- 
tities d and oj. They provide detailed information 
about the accuracy of structure-optimization and frozen- 
phonon calculations performed in this paper. 



VI- 



the two alternatives 



6. Vacuum Size 



The first set of test calculations serves to determine 
the appropriate size of the vacuum region. The results 
presented in Table ^ show that an equivalent of four 
substrate layers (= 8.8 A) is sufficient to decouple the two 
substrate surfaces accross the vacuum. Both relaxation 
parameter d and the frequency lo vary only within the 
given error range upon further increasing the vacuum 
region. 



c. LAPW Parameters 



Another important parameter is the wavefunction en- 
ergy cutoff 5^(^"^)^- It has a decisive influence on the 
computer time. With res pect to the structure parame- 
ter d our results in Table 



possible value is 2^{K' 



wf 



VI1| indicate that the smallest 
= 12 Ry. However, for hu) 
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TABLE VIII. Mo(llO) double-layer: Variation the poten- 
tial paremters ^(G^"*)^ and The notation is equivalent 
to Table [v^. 

^{G^°'f d [%do] huj [meV] 



[Ry] 


E 


F 


E 


F 


64 


-2.0 


-2.3 


37.3 


35.6 


81 


-1.9 


-2.3 


36.6 


35.5 


^ 100 


-2.0 


-2.4 


37.3 


35.6 


121 


-1.9 


-2.3 


37.0 


35.6 


144 


-1.9 


-2.3 


37.0 


35.6 



P°* d [%do] Tiu [meV] 

E F E F 

2 =19 =T?7 37!5 346 

^ 3 -2.0 -2.4 37.3 35.6 

4 -2.0 -2.5 37.1 35.2 

5 -2.1 -2.6 37.1 35.3 



already a calculation with ^{K^^^Y = 10 Ry provides 
converged results. 

The two potential parameters and are 

of little influence on t he calculated quantities. This can 
be seen in Table VIII. The values of the relaxation pa- 
rameter d and the frequency lo remain relatively sta- 
ble within the parameter range 64 Ry < < 
144 Ry. Both quantities are more sensitive to changes in 
the {Urn) expansion of the MT potential. However, a 
value of /P°* = 3 seems to be sufficient. 

In general we note that energy (E) and force (F) cal- 
culations demonstrate a similar convergence behavior. 
However, the magnitude of the relaxation in the force cal- 
culations is overestimated by about Ac? — 0.?)%<F°™ « 
0.01 A; for the frequencies one finds a deviation of about 
fiAuj = h%fujj'^°^'^ K, 1.8 meV from the converged total 
energy values. Those differences are probably due to nu- 
merical inaccuracies in the force calculation. 



d. kii Point Set 



Similar calculations for a W(llO) double-layer system 
lead to a converged LAPW parameter set which is com- 
parable to the one obtained for Mo. Only the MT po- 
tential parameter Zp°' has to be chosen higher: Zp°' = 4. 
slightly different. The agreement between force and to- 
tal energy calculations (Ad = « 0.02 A and 
Ao; = b%uj'^°™) is nearly as good as in the case of Mo. 
The following calculation is aimed to determine a two- 
dimensional k|| -point set which can be used for the calcu- 
lation of the atomic and electronic structure of Mo (110) 
and W(llO) surfaces (see Table |^). Because of the size 
of (110) surface slab systems it is important to keep the 
number of k|| -points as small as possible. For the cal- 
culations we use the LAPW parameter set given in Ta- 
ble I and vary the k|| -point set for either the valence or 
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TABLE IX. W(llO) double- layer: Accuracy of the calcu- 
lation in dependence of the number of ky -points for the (a) 
valence electrons and (b) semicore electrons. The respective 
fixed k||-point set is made uo by one special ky-point in the 
irreducible part of the SBal3. The arrows indicate the con- 
verged k||-point set used for the ( 110) surface calculations. 
The notation is equivalent to Table VI. 



(a) valence electrons 



# d [Vodo] hoj [meV] 



k|| -points 


E 


F 


E 


F 


36 


-8.5 


-9.5 


20.5 


21.1 


49 


-9.0 


-10.1 


23.8 


23.1 


^ 64 


-8.4 


-9.3 


19.2 


20.3 


81 


-8.4 


-9.4 


18.9 


20.4 



(b) semicore electrons 



# 

k|| -points 


d [%do] 




hu} [meV 




E 


F 


E 


F 


4 


-1.0 


-1.8 


27.4 


26.2 


^ 9 


-1.0 


-1.9 


27.4 


26.2 


16 


-1.0 


-1.9 


27.3 


26.2 


25 


-1.0 


-1.9 


27.4 


26.2 


36 


-1.0 


-1.9 


27.4 


26.2 



the semicore electrons. These sets always consist of uni- 
form two-dimensional point meshes. Converged results 
are obtained by using 64 k|| -points for the valence and 9 
k|| -points for the semicore electrons. 



3. Phonons 



a. LAPW parameter 



Our final series of tests is meant to determine the ac- 
curacy of the frozen-phonon calculation. In order to keep 
the computer time for those elaborate calculations within 
pleasant limits one has to keep the number of ky -points 
and basis functions as small as possible. As a sample case 
we study the S point Rayleigh phonon of the H/Mo (110) 
system. 

The calculations are performed as follows: We dis- 
tort the atoms of the first substrate layer by tiz{S) = 
0.08 A according to the pattern in Fig. ^ and deter- 
mine the total energy and force changes AiJ*°* and AFi^. 
These quantities have the same convergence behavior as 
the phonon frequencies: Therefore, if AE^°^ and AFi^ 
are stable with respect to a particular LAPW parameter 
so are the phonon frequencies. Converged results for Mo 
as well as for W are obtained by using the following set- 
ting: wavefunction plane- wave cutoff |^(i4r'^^)^ = 10 Ry; 
pot _ 2 .j-jjg j-^^ j^-j expansion of the MT potential; uni- 
form k|| -point mesh consisting of 56 points; five substrate 
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TABLE X. Frozen-phonon distortion for H/Mo(110): 
Change of tiie total energy Ai5'°* and the force AFiz ver- 
sus the number of ky-points. All ky-point sets are uniform 
meshes which include the symmetry point S' of the backfolded 
SBZ (see Fig. |ll|) . Given are also the calculated values for the 
clean Mo(llO) surface (with 56 k||-points). 



# k||-points 


A£'°' [meV] 


AFi, [meV/A] 


16 


110 


672 


56 


136 


756 


80 


144 


797 


120 


135 


790 


232 


136 


777 


clean 


161 


938 



layers; a vacuum region equivalent to four substrate lay- 
ers; XC potential from Refs. 
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b. k|| Point Mesh 

In order to accurately describe the electronic structure 
and hence the electron-phonon coupling for the S point 
phonon one has to provide a detailed k||-point set. For 
instance, for the study of the (2x1) reconstructions of 
the diamond (111) surface Vanderbilt and Louie employ 
a k|| -point mesh which becomes logarithmically denser 
close to the zone boundarytj. One cannot rely on the 
results for the two M o (110) and W(llO) double-layer 



slabs in Section A 2 d because in these systems nesting 
effects are totally unimportant. 

We perform tests for the H/Mo(110) system using up 
to 232 k|| -points in the SBZ and find that a uniform k||- 
point mesh of 56 points is sufficient in order to obtain con- 
verged results (see Table ^). For all those frozen- piionon 
calculations the Monkhorst-Pack k|| -point meshesc3 are 

shifted in order to include the S point of the (1x1) unit 



cell (see Fig. 11). For the k||-point summation we em- 
ploy a Fermi smearing of fcBToi ~ 68meV. Only for the 
232 k||-point calculation listed in Table ^ a smaller value 
of k^Tci ~ 14meV is used. However, this parameter is 
apparently uncritical because the energies are obtained 
for the /cB^ei ^ K limit as discussed in ReL E2. 



c. Basis functions 

Next, we check whether the results are sensitive to an 
increase in the wavefunction cutoff, i.e., whether the cho- 
sen value iL_(j;^wf^2 _ 10 Ry is large enough. The data 



listed in Table XI demonstrates that this is indeed the 
case. The results remain stable for the higher cutoff 
|i(if-f)2 = 12Ry. 
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TABLE XI. Frozen-phonon distortion for H/Mo(110): 
Change of tiie total energy A£'°' and the force AFi^ with 
respect to the wavefunction plane-wave cutoff A_(ifwf^2^ 

^{K-^y [Ry] AE'°' [meV] AFi, [meV/A] 

10 136 756 

12 137 767 



TABLE XII. 


Frozen-phonon disto* 


tion foi^/Mo(110): 


Variation of the 


XC potential: LDAEIE 


i ggaeI. 






XC Potential 


A£*°* [mcV] 


AFi, [meV/A] 


LDA 


136 


756 


GGA 


141 


764 



d. LDA versus GGA 



For some systems involving hydrogen the results of 
a total energy calculation are not reliable if the LDA 
is used. In those cases the bonding — especially be- 
tween hydrogen atoms — is only described accurately if 
gradient corrections are consideredO. In Table KIl we 
compare results obtained with the generalized gradient 
approximation (GGA)cJ'C3 and those found within the 
LDA. The latter is identical with the 56 k|| -point calcu- 
lation in Table-^. There, the XC potential of Ceper|i 
ley and Aldeic^parametrized by Perdew and Zunge: 
is used. The small energy and force differences of 5 meV 
and 8meV/A, respectively, clearly shows that our adsor- 
bate systems are rather insensible to the XC formalism 
chosen. 
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